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The evolution equations for the generalized microscopic phase densities are introduced. The 
evolution equations of average values of microscopic phase densities are derived and a solution of 
the initial-value problem of the obtained hydrodynamic type hierarchy is constructed. 



For the kinetic description of plasmas the method of 
microscopic phase density is successfully used P, 0, 
The aim of the present paper is to derive the evolution 
equations for the generalized microscopic phase densities 
and their average values, i.e. the hydrodynamic type 
equations, in a more consistent way, and construct solu- 
tions of the corresponding initial- value problems. 

We consider the system of a non-fixed (i.e. arbitrary 
but finite) number of identical particles with unit mass 
m = 1 in the space R 3 (nonequilibrium grand canon- 
ical ensemble). Every particle is characterized by the 
phase space coordinates Xi = (qt,Pi), i.e. by a posi- 
tion in the space qt £ R 3 and a momentum pi 6 R 3 . 
A description of many-particle systems is formulated in 
terms of two sets of objects: by the sequences of observ- 
ables A = (A ,Ax(x 1 ),...,A n (xi,...,x n ),...) and by 
the sequences of states D — (1, Di(xi), . . . , D n (x\, . . . , 
x„), . . .). The average values of observables determine 
a duality between observables and states. As a conse- 
quence, there exist two approaches to the description of 
the many-particle system evolution, namely those con- 
cerning the evolution of observables or the evolution of 
states: 
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where (l,£)(0)) = J2n=oi\ I dx i ■ ■ -dx n D n (0) is a nor- 
malizing factor (grand canonical partition function) . The 
sequence D(t) = (l,Di(t,xi), . . . ,D n (t,xi, . . . , x n ), . . .) 
of probability densities of the distribution functions 
D n (t) is a solution of the initial- value problem of the 
Liouville equation. The sequence of observables A(t) = 
(Aq, Ai(t, xi), . . . , A n (t, xi, . . . , x n ), . . .) is a solution of 
the initial- value problem of the Liouville equation for ob- 
servables. If A(Q) is the sequence of continuous functions 
and D(0) is the sequence of integrable functions, then 
functional (JTJ) exists. 

An equivalent approach of the description of evolution 
of many-particle systems, that enables to describe 
systems in the thermodynamic limit, is given by the 
sequences of s-particle (marginal) distribution func- 



tions F(t) = (l,F 1 (t,xx),...,F a (t,Xi,...,x a ),...) 
and s-particle (marginal) observables G(t) = 
(G ,Gi(t, xx),..., G a (t,xi,...,x s ),...). The se- 
quence F(t) is a solution of the initial- value problem of 
the BBGKY hierarchy and G(t) is a solution of 

the initial-value problem of the dual BBGKY hierarchy 
0. In that case, the average values of observables at 
time moment tel are determined by the functional 

(A)(t) = jr i - f dx 1 ...dx a G 8 (0)F s (t)= (2) 
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Thus, the sequence of marginal observables G(t) in terms 
of the sequence A(t) is defined by the formula 
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where Y = (xi, . . . , x s ), s > 1, and the sequence F(t) of 
marginal distribution functions is defined in terms of the 
sequence D(t) as 



F a (t,xi, ...,x a ) = 

00 1 . 
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We remark that, in the case of a system with a fixed 
number N of particles (nonequilibrium canonical en- 
semble) observables and states are the one-component 
sequences, respectively, A^ N ^ = (0, . . . , 0, An, 0, . . .), 
flW = (0, . . . , 0, D N , 0, . . .). Therefore, the formula for 
average value ([1]) reduces to the expression 

(AW) = (l, D^y 1 fdxx... dx N A N D N , 

where (l, #>) = J dx x . . . dx N D N is a normalizing fac- 
tor (canonical partition function). 

We introduce the observables known as the 
microscopic phase densities of the system of a 
non- fixed number of identical particles. Let 
N(t) = (NW(t),...,NW(t),...), where N^(t) = 
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(Q,...,0,Njf\t),...,N} k >{t),...), k > 1, is the 
sequence of microscopic phase densities of fc-ary type 

(f) = jv« (t, a*, ...,*„) = (5) 

n A; 

where $ is the Dirac <5- function, £i , . . . , £fc are the macro- 
scopic variables £j = (i>i, 7"j) £l 3 xll The set of func- 
tions {Xi(t, xi, . . . , x n Y\ ._ 15 n > k > 1, is a solution of 
the Cauchy problem of the Hamilton equations for n par- 
ticles with the initial data x\ , . . . , x n and with the Hamil- 
n 2 » 

tonian iJ„ = J] ^ + 2 — where - g^) 

i=l i<j = l 

is a two-body interaction potential. 

For example, if /c = 1, i.e. in the case of an additive- 
type observable, we have the microscopic phase density 

n 
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Microscopic phase densities ([5]) are the solutions of a 
sequence of the Cauchy problems of the Liouville equa- 
tions for observables 
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with the initial data (1 < fc < n) 
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where the brackets (•, •) denote a scalar product of vec- 
tors. 

We note that solution |[SJ of Cauchy problem ©-(O 
defines the one-parametric group of operators R 1 3 i n 
S n (t)N n (0), i.e. 

N<p(t, & U--r : x n ) = S„(i)/V«(0), (8) 

where N^ k) (0) is the microscopic phase density 0. 

In terms of variables £i , . . . , , the sequence of Liou- 
ville equations © for microscopic phase densities ([5]) is 
represented as the BBGKY equations set with respect 
to the arity index k > 1, while it is a sequence of equa- 
tions with respect to the index of the number of particles 
n > k. Indeed, we have 

i—1 

+ E(|:^-^|:))^) + 

+ E/^x<^*(r i -r M . 1 ) > A>^)(*). 

i— 1 



If fc = n, it is the Liouville equation. Using such a rep- 
resentation of ©, we can directly derive the evolution 
equations for average values |lj of microscopic phase 
densities ([5]). We observe that in the case of a system 
of a fixed number TV of particles, equations ([9]) are the 
hierarchy-type equations set with respect to the index 
k > 1 , and the Liouville equation with respect to a num- 
ber of particles. 

The fact that the microscopic phase densities ([5]) are 
exactly governed by the BBGKY-type hierarchy of equa- 
tions © is closely connected with the structure of the hi- 
erarchy. Indeed, it is known [6( that the BBGKY hierar- 
chy for the states has the explicit solution - the marginal 
pure state which is the sequence of functions of the phase 
space variables similar to microscopic phase densities ([5]) 
as functions with respect to the " macroscopic variables" . 

We introduce the sequence of marginal observables 
G(t) = (G« (*),..., (<),...) of /c-ary type (t) = 

(0, . . . , 0, G [k) (t), . . . , G { s k) (t),...) defined by © through 
microscopic phase densities ©. For example, according 
to ([3]) at the initial time moment t = the sequence of 
marginal additive-type microscopic phase densities has 
the form £?W(0) = (0, . . . , d^-x^), 0, . . . ). Correspond- 
ingly, the sequence of marginal observables of fc-ary type 
microscopic phase densities ([5]) is given as follows: 

G«(0)= (10) 
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= (0,...,0, Y, II*(6-a5i,),0,...)- 

il#...^ifc = l 1=1 
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Let Y = (xi,...,x s ) and (xi, . . . , V, . . . , x s ) = 
(xi,...,Xj-i,Xj+i,...,x 3 ) = Y \ Xj. The 
marginal microscopic phase densities G^ (t) = 

Gi fc) (i,^i, . . . ,£k',xi, . . . ,x s ) of every fc-ary type are 
governed by the initial-value problem of the dual 
BBGKY hierarchy 

!g?>(,)-(j>.£>- an 
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with the initial data 

G^(i) | t=0 =G«(0), s >fc>l. (12) 

As a case in point, we adduce the first equation of hier- 
archy (fTTj) 

i—1 
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For the marginal additive-type microscopic phase density, 
the first two equations have the form 

2 
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To construct a solution of the dual BBGKY hierarchy 
(|11[) . we introduce some preliminaries. On continuous 
functions we introduce the nth-order cumulant, n > 1, of 
the groups of operators 

% n (t)=% l {t,X)= (13) 

= e (-i) ipm (ipi-i)! n s \xM 
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where is the sum over all possible partitions P of 
the set X = (xx, . . . ,x n ) into |P| nonempty mutually 
disjoint subsets Xi C X, the operator S|x 4 |(i) is defined 
by formula ©. 

On a continuously differcntiable function g\ — gi(xi) 
the generator of the first-order cumulant is defined by the 
operator 
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lim-(2li(i,a;i) - l) gi = (p t , ^—)gi- 
t—>o t v ' aqi 

In the case n = 2, we have 

1 & d 

Yua- t Mt,x u x 2 )g 2 = - ]T ( — - — ><*. 



If n > 2, as a consequence of the fact that we consider 
a system of particles interacting by a two-body potential 
the limit 

lim j% n (t)g n = 0. 

holds. 

Let Y = (xi, . . . , a; s ), X = Y \ {xj 1 , . . . , Xj a _ n }. For 
continuous functions in the capacity of initial data, a 
solution of Cauchy problem (fTTj) - (fT2|) is an expansion 
over particle clusters whose evolution are governed by 
the corresponding-order cumulant (semiinvariant) of the 
evolution operators of finitely many particles 

Gi k \t,Y) = j2^ E W*. ( 14 ) 

n=0 ' j 1 =i...=i 3n =l 

(Y \ A-)GW n (0, y \ {a* , . . • , x jn }), 



where the evolution operator 

% 1+n (t,(Y\X) 1 ,X)= ]T (-!)' Phl 

P:{(n^)i.^}=U^i 

x(|P|-l)! ft S lXil (t,Xi) 

XiCP 

is the (1 + n)th-order cumulant [8[ of groups StxMt) of 
operators |8]l. ^ p is the sum over all possible partitions 
P of the set \ (Y \ X)i,X} into |P| nonempty mutually 
disjoint subsets X, C {(Y \ X)) V X}. The set (Y \ X)x 
consists of one element of Y\X, i.e. the set Y\X = 
{xj 1 , . . . , Xj ■ _ } is a connected subset of the partition P 
(|P| = !)■ 

For the additive- type microscopic phase density (|10[) . 
we derive 

GW(t,Y)=* a (t i Y)^26(Z 1 -x j ). (15) 
i=i 

Then in terms of variables £i , . . . , £/. , the first equation 
of hierarchy (jlip for the additive-type microscopic phase 
density (fT5"|) takes the form 

^(Ui;^,..,^) = 

d 

= -(v 1 ,—)GW{t,£i;x 1 ,...,x s ) + 
+ y rf6(^1>(ri-r 2 ),— )Gi 2 )(t,6,6;a ; i,.--^ s ), 
with the initial data 

s 

G^sHt, ^Xi, ...,X S ) | t=0 = E^^ 1 ~ X i) S s,l, 

i=l 

where 5 Si i is the Kronecker symbol, s > 1. 

In a similar manner for the marginal microscopic phase 
densities of fc-ary type G«(t) = (0, . . . , 0, G^\t), 
G s (t), . . . ) we derive 



Of 



(16) 



+ £ / da . +1 (|-<i>(r l -r fc+1 ),^-}Gf+ 1 )(t) 

-i— 1 

with the initial data 

s k 

Gi k Ht)\t=o= E ( 17 ) 

Here, 1 < r < s, and if fc = s, the marginal microscopic 
phase density G^(t) is governed by the Liouville equa- 
tion. 
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Thus, in terms of variables £1, • • • , £fc the dual BBGKY 
hierarchy (|16p for marginal microscopic phase densities 
(fl"4|) is represented as the Bogolyubov equations set with 
respect the arity index k > 1, while evolution equations 
(fl~6)) have a structure of a sequence of equations with 
respect to the index of the number of particles s > k. 

We introduce the evolution equations for average val- 
ues of marginal microscopic phase densities (1141) . For the 
k-aiy type microscopic phase density, according to ((2|) 
and lfl"j]). from the dual BBGKY hierarchy (JTSJ) we de- 
rive the following hydrodynamic type equations for their 
average values © : 

!<*»>(«)= (-£<«■.£>+ <i8) 

2 — 1 

+ E<^*(ri-ri) 1 ^»(G«>(t) + 

+ E / ^ + i(^$(n-r fc+1 ),^)(G^)}(0, 

2 — 1 * 

with the initial data 

<GW)(t,£ 1 ,...,a)| f=o HG (fc) >(0), *> 1.(19) 

Due to functional ([2J, initial data (fT9)) are given as the 
functions (G^)(0, £ 1} . . . , &) = F fe (0, 6, • ■ ■ , &)> where 
^(0, £i, . . . , £ s ) is the value of initial marginal state at a 
point £i, . . . ,£ s . The hierarchy of equations (18) with re- 
spect to the arity index we refer to as hydrodynamic type 
equations since it describes the evolution of the average 
values of the microscopic phase densities. 

We note that, according to the definition of functional 
© and ©, the equality (G (fc) )(t) = (N^)(t) holds in 



the case of finitely many particles. In the thermodynamic 
limit, the value (N^)(t) tends to (G^)(t), i.e. to the 
solution of Cauchy problem p^ |) - ([T9|) . 

A solution of Cauchy problem (fT8|) - (ri9|) is defined 
by the expansion over the arity index of the micro- 
scopic phase density, whose evolution is governed by the 
corresponding-order cumulant of the evolution operators 
similar to (|8|), namely 

oo „ 

<<?<*>>(*,&,...,&) = £--1 <% k+1 ... (20) 

n=0 ' J 

■ ■ ■ ^+„2ti+„(-i, n.St+i, ■ ■ ■ xo), 

where (G ( - k+n ^)(0) are initial data (JTTJ) and 
2ti+n(-t, Yi,£ k +i, ■ ■ ■ ,£k+n) is the (1 + n)th-oidei 
cumulant of groups of evolution operators similar to JHJ). 
For integrable functions (G^ fe ))(0) series (|2T))) converges 
for small densities. We note that, if we apply the 
Duhamel formula to cumulants of groups of operators 
similar to ([5]), solution expansion (|20[) reduces to the 
iteration series of the hydrodynamic- type hierarchy (|18[) . 

In summary, the evolution of the generalized marginal 
microscopic phase densities (fit))) is described by the 
initial- value problem of the dual BBGKY hierarchy (jlip . 
Their average values © are governed by the initial- 
value problem of hierarchy (|18[) of the hydrodynamic- 
type equations which has the structure of the BBGKY 
hierarchy with respect to the arity index of the gener- 
alized microscopic phase densities. The solutions of the 
Cauchy problems of such hierarchies (fTT|) and (fl~8|) are 
represented by the expansions (TT4")) and (|20p correspond- 
ingly. 
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